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Abstract 

We present a pedagogical case study how to combine micro-causality and uni- 
tarity based on a perturbative approach. The method we advocate constructs an 
analytic extrapolation of partial-wave scattering amplitudes that is constrained by 
the unitarity condition. Suitably constructed conformal mappings help to arrive 
at a systematic approximation of the scattering amplitude. The technique is illus- 
trated at hand of a Yukawa interaction. The typical case of a superposition of strong 
short-range and weak long-range forces is investigated. 



1 Introduction 

Recently a novel scheme for studying hadronic interactions beyond the thresh- 
old region was introduced by two of the authors pp. The main objective of 
that work is a controlled realization of the causality and unitarity condition 
in a perturbative application of the chiral Lagrangian. The starting point are 
partial-wave dispersion relations. A generalized potential is constructed from 
the chiral Lagrangian in the subthreshold region and analytically extrapo- 
lated to higher energies. The partial-wave scattering amplitudes are obtained 
as solutions of non-linear integral equations. 

The purpose of the present letter is an illustration of the method for a schematic 
system where the exact solution is known. This will shed further light on 
its usefulness. We consider non-relativistic Yukawa interactions of various 
strengths and ranges. Though a quantitative description of for instance the 
pion-nucleon or nucleon-nucleon system is not justified in terms of static 
Yukawa potentials there are essential features grasped by a Yukawa ansatz. 
The long-range part of the nucleon-nucleon force can be understood in terms 
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Fig. 1. Convergence domain 
of the Born series for the 
s-wave Lippmann-Schwinger 
equation ([!]) with Q as 
a function of the coupling 

■ strength g/g c - The dashed 

curve locates the branch point 
"^■^^ caused by the 2nd order Born 

term. 



of a static Yukawa potential derived from the pion-exchange process. The im- 
portant feature is a left-hand cut in the partial-wave scattering amplitudes 
very close to the physical region. In this respect the pion-nucleon system is 
similar, as its long-range force generated by the u-channel nucleon exchange 
process defines close-by left hand cuts as well. 



2 Expansions with a Yukawa potential 



We consider the quantum mechanics of two particles with masses mi and 
m 2 denned by a single Yukawa potential. The system is characterized by the 
reduced mass parameter m = mi m 2 /(mi + m?)- The scattering process in a 
given partial wave with angular momentum I is determined by the t-matrix, 
satisfying the nonrelativistic partial- wave Lippmann-Schwinger equation [2]: 

<* i ^ , i*> = m ,*> + ^ s: *- «• ( ^™f ifc> , (i) 

with q, k, k' being the center of mass on-shell and initial and final off-shell mo- 
menta respectively. The Yukawa potential projected onto angular momentum 
reads 



,,„„,,, g „ / k 2 + k' 2 + fi 2 \ ^ . , 1, x + 1 

{kKm= 2*k Q <{ 2Vk - )< *W = -log — , 

Q,(x) = if,(x) log|±I-£ [p.-iWPhW, (2) 



k=l 



with the conventional Legendre polynomials Pi(x). The Yukawa mass /x charac- 
terizes the range of the interaction. For simplicity and definiteness we assume 
nonrelativistic kinematics. 



2 



The integral equation ([T]) can be expanded in powers of the coupling constant 
g, defining the Born series. The expansion converges for certain values of g 
and q 2 [HE]- The area of convergence in the plane of q 2 /fj, 2 and \g/g c \ is shown 
in Fig. 1 for the s-wave potential. For coupling constants, \g\ < g c , smaller 
than the critical one 



C ~1.68-^-, (3) 
2 m 

the whole region — oo < q 2 < +oo is perturbative. Within the shaded area 
the Born series does not converge. If the coupling constant exceeds its critical 
value the scattering amplitude is perturbative only either for sufficiently small 
or large q 2 as shown in Fig. [2} In the case of attraction with g < the first 
bound state is formed at g = — g c [5J. An important observation is that even 
if perturbation theory fails to converge close to threshold, it does converge for 
sufficiently large negative q 2 . 

We confirm the convergence properties by an explicit computation of the Born 
series to third order. In Fig. [2] the scaled on-shell amplitude 



U(z) =2mfi(k'\ U(q 2 )\k) 



q 2 

with z = \, (4) 

--k=q 



is shown for the particular choice \g\ = 3g c /2 in an s-wave with I = 0. The 
dimensionless amplitude ii(z) is a function of the ratios z = q 2 /fi 2 and g/g c 
only. While the left-hand panel in Fig. |2]follows with positive and repulsive, the 
right-hand panel with negative and attractive coupling constant. The figure 
illustrates convergence for large \q 2 \. In contrast the near-threshold region 
cannot be described by the Born series. Most striking is the bound state 
structure which arises for negative coupling constant. Higher partial waves 
can be treated similarly to the s-wave, but the speed of convergence for them 
is in general higher, the critical coupling constants are significantly larger then 
the one for the s-wave. 

How to construct a systematic approximation of the partial-wave scattering 
amplitude for the non-perturbative case with \g\ > g c 7 While for the given 
Yukawa potential this is a rather academic question this is not the case for 
realistic systems that can typically not be solved exactly. The idea is to con- 
struct an expansion that preserves the causality and unitarity constraints of 
the scattering amplitude pp. 

For our further consideration it is useful to introduce yet another dimensionless 
scattering amplitude. While the amplitude Q illustrates the scale dependen- 
cies of the system, it suffers from kinematical constraints at threshold. The 
angular momentum barrier leads to 
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Fig. 2. S-wave scattering amplitude for \g\ 
for z ~ . 



q 2 V 



3_3c 
2 



We introduce an on-shell partial-wave scattering amplitude 



(5) 



TM 2 



2 m J A 2 + q 2 



'A 2 + q 2 Y 



{k'\tM 2 )\k) 



k'=k=q 



(6) 



with the scale parameter A. By construction the dimensionless partial- wave 
amplitudes are finite at threshold with q 2 = 0. Moreover, for any finite value 
of A the amplitudes vanish with large momentum q 2 as 



TM 2 



m q 



2v? 



for 



g 2 »A 2 . 



The unitarity constraint takes the simple form 



(7) 



ZT l (q 2 ) = \T l (q 2 )\ 2 p l (q 



pM 2 



A 2 + q 2 



(8) 



illustrating that the parameter A characterizes the scale at which the phase- 
space function pi(q 2 ) approaches its asymptotic value 1. If the unitarity condi- 
tion (JsT) is satisfied the partial-wave scattering amplitude can be parameterized 
by the phase shift Si(q 2 ) with 



Ti{q 2 ) - l.( ■ I 



(9) 



The analytic properties of the scattering amplitude generated by a Yukawa 



4 



potential is well understood (see e.g. [6]). A key observation is the known 
analytic structure of the generalized potential U[(q 2 ) introduced by 



UM 2 ) = TM 2 ) - f ^ ,/ (g ? . I W 2 )l 2 , (io) 

Jo 7r q' 2 + \r M q' 2 — q 2 — te 

where we insist on a matching scheme (see e.g. [TflT] ). At the matching scale 
q 2 = —[1m the generalized potential and scattering amplitude take identical 
values. While the scattering amplitudes Ti(q 2 ) possess a branch point at q 2 = 0, 
reflecting the unitarity cut, the generalized potential is analytic at q 2 > 0. The 
potential is characterized by left-hand cuts only. For a Yukawa potential it has 
branch points at q 2 = —A 2 and q 2 = — \ (n p) 2 with n = 1, 2, ... [BJ. 

A crucial step is to apply a perturbative expansion in powers of the coupling 
constant g not to the scattering amplitude but rather to the generalized po- 
tential Ui(q 2 ). Such an expansion converges even for \g\ > g c if the matching 
point \im is chosen outside the nonperturbative domain of Fig. [2j Previous 
but somewhat different studies along this line are [8f9fTU] . 

We focus on s-wave scattering with \g\ — | g c and choose a matching scale 
H 2 M = 10 fi 2 , that lies well below the boundary of the nonperturbative region. 
The parameter A will be varied. For a given truncation of Ui(q 2 ) we use (10) 



to reconstruct the corresponding approximate scattering amplitude. The non- 
linear integral equation can be solved by means of the N/ D technique [TT|T2] . 
The scattering amplitude is decomposed into the form 



T(g 2 ) = ^ (11) 

where the function N(q 2 ) has only left-hand singularities, and D(q 2 ) has only 
right-hand singularities. These functions obey the following system of linear 
equations 



N(q 2 ) = U(q 2 ) + 



-/2 



g 2 + / 4 N(q' 2 )p(q' 2 )(U(q' 2 )-U(q 2 )) 



o n q' 2 + Hm q' 2 - q 2 



fl(rt=1 _ / v ^mm . (12) 

' 7r q' 2 + n 2 M q' 2 — q 2 — te 



In Fig. |3]we present the result of our numerical simulation of (12) for various 
truncations of the generalized potential. We verified that indeed for all consid- 
ered cases a solution of the non-linear integral equation (10) was generated. In 
order to offer a more quantitative analysis we provide the scattering phase shift 
rather than the s-wave scattering amplitude. The figure demonstrates that the 
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Fig. 3. Scattering amplitude for g = | g c (left) and g = —\g c (right), 
bands are implied by a variation of parameter A with 3 [i < A < 9 \i. 



The various 



perturbative evaluation of the generalized potential recovers the exact phase 
shift with increasing precision uniformly in energy. This is in contrast with the 
Born series, which does not recover the scattering amplitude in the threshold 
region. Already a truncation of the generalized potential accurate to cubic or- 
der reproduces the exact phase shift quite accurately. The result is insensitive 
on the parameter A. The spread implied by a variation 3/i<A<9/iis indi- 
cated by the bands. As expected the band width is decreasing systematically 
as the order of the truncation increases. 



For the attractive case there is a bound state in the system. In the solution of 



(12) it appears as a zero in the D-function and therefore the scattering ampli- 
tude develops a pole at q 2 = —q 2 B - The binding energy is E B = q 2 B /{2m). In 
this case the solution of (12) does not constitute a solution of the nonlinear 



integral equation (10). The truncated potential does not have a pole contri- 



bution. However, the solution of (12) is the solution of a related system where 
the pole term 



{q 2 + Ai)gl ( 13 ) 



(q 2 + 2mE B ) (2mE B -fi 2 M ) 
is added to the approximated potential. It is straightforward to show that the 



set of corresponding N/D equations (12) defines a scattering amplitude, which 



coincides with the one in the absence of (13) provided that the pole residuum 
g B is dialed properly. 

Though the result of Fig. [3] is encouraging it does not yet permit a direct 
generalization to realistic systems. In effective field theories, such as Chiral 
Perturbation Theory one can perform a perturbative expansion only for small 
momenta. At first it appears impossible to arrive at a generalized potential 
that can be trusted in the region needed to reconstruct the partial-wave scat- 
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tering amplitude via the non-linear integral equation (10). The crucial obser- 
vation is that the potential is required for q 2 > only. This together with the 
absence of branch points at q 2 > in the potential permits the set up of a 
controlled approximation strategy for the generalized potential. 

We split the potential into an inside and outside part 

U(q 2 ) = U inside (q 2 ) + U outsidc (q 2 ) , (14) 

where [/inside (? 2 ) contains the cuts, closest to the physical region and therefore 
have the largest effect on the energy dependence of the generalized poten- 
tial. Moreover, the inside part of the potential can be obtained with ease 
in particular within effective field theories. The strength of an effective field 
theory is its capability to compute threshold properties quite accurately and 
systematically. All contributions from the more distant cuts are subsumed 
in £/ utside(g 2 )- Typically, the more distant cut structure are quite difficult to 
compute. In effective field theories most frequently this is impossible. 

In the case considered it is natural to include into Ui nsi( ie(q 2 ) the cut from 
— fi 2 < q 2 < — ^ that is determined exclusively by the first Born term. It 
takes the form 



U insidc (q)=gm J ^— ^ 4 ^ /2 . (15) 
-p? 

In turn the outside part of the potential is characterized by branch cuts in the 
interval (— oo, — /j, 2 ). 

If the Taylor expansion of £/ utside(<? 2 ) around q 2 = is known, we can use 
the method of conformal mapping to analytically extrapolate the generalized 
potential onto the whole positive real axis with q 2 6 (0, oo). Even though in our 
case the Taylor series converges for \q 2 \ < /i 2 only, the knowledge of the Taylor 
coefficients permits an accurate reconstruction of the outside potential in the 



region required in (12). This is achieved in application of the transformation 



2 fi- yjq 2 + fi 2 

) = , / 2 i 2 ' ( 16 ) 

/i + v <r + fi z 

which maps the cut complex plane onto the unit circle around £(0) = 0. The 
generalized potential can be reconstructed with 



U(q 2 ) = U inside (q 2 ) + J2C k [aq 2 )} • (17) 

A:=0 
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Fig. 4. Scattering amplitude for g = | g c (left) and g = —\g c (right). The various 
bands are implied by a variation of parameter A with 3 [i < A < 9 \i. 

If the first n coefficients in a Taylor expansion of outside (9 2 ) are known, one 



can compute Cy. for k < n. By construction the expansion (17) converges in 
the cut complex plane. 



We studied the convergence properties of (17) for the examples considered 
above. An expansion to the exact potential of the form (17) was performed. 
The results are collected in Fig. |4| A rapid convergence is observed for the 
repulsive and attractive system. An accurate reproduction of the exact phase 



shift is achieved considering the first three or four terms in (17) only. Like in 
Fig. [3] the sensitivity of our result on the parameter A is illustrated by bands. 
The band width decrease systematically with increasing order of the expan- 
sion. For the sake of clarity only two approximative results are confronted in 
Fig. [4] with the exact solid line. 



With (17) we present an alternative expansion scheme for the generalized po- 
tential that appears more practical in realistic applications. The coefficients Ck 
may be computed approximatively in perturbation theory or may be adjusted 
to empirical data directly. 



3 Expansion after renormalization 



The system described in the previous section shows some partial analogy with 
effective field theories which may be constructed to be perturbative within the 
Mandelstam triangle. In contrast to the Yukawa toy model in a realistic sys- 
tem the matching scale hm cannot be chosen arbitrarily low. Nonperturbative 
effects from the u- or t-channel restrict the perturbative domain from below. 
The matching should be performed within the Mandelstam triangle. In order 
to account for this requirement we set Hm — A* in what follows. This particular 
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choice implies a matching right at the t-channel unitartity branch point. 

A further important issue is the intricate interplay of short and long-range 
forces. Typically the long-range forces are rather week and perturbative but 
the short-range forces are strong and nonperturbative. For example in nucleon- 
nucleon scattering the central part of the pion-exchange potential has a cou- 
pling constant that is three times smaller than its critical value [5]. On the 
other hand there is often a strong short-range contribution, which can be 
responsible for the formation of bound states or resonances. The nature of 
short-range interactions is not resolved by effective field theories. They get 
renormalized as to meet the physical conditions. We model this situation by 
considering an s-wave system defined by a superposition of two Yukawa po- 
tentials 



{k'\V \k) = (k'\V L \k) + (k'\V s \k) 



g L (k + k') 2 + fi 2 L g s (k + k') 2 + fi 2 s 

i0 S 7.A9 , ..2 + i 10 g 



4k'k ~° (k - k') 2 + & 4k'k {k — k') 2 + /4' 
where the second term is short range in the sense that fig ^> fiL- 

We renormalize the short-range force as to arrive at a closer correspondence 
to effective field theories. In an initial step we treat the short-range part of 
the potential exactly but seek a perturbative treatment of the long-range part. 
For this purpose it is useful to introduce a compact matrix notation for ([TJ. 
We write 



1-(V l + V s )G] l [V L + V £ 

oo 

t s + {l+t s G)V L J2[(G + Gt s G)V L 



n=0 



{l + Gt s ), 



(19) 



where tg = Vg + Vs G t$ is the full t-matrix in the presence of the short-range 
potential alone. The crucial observation is that a contribution which involves 
n long-range potentials Vl and k short-range structures t s is characterized by 
its left-hand cut starting at 



{nfi L + kfi s ) 2 ■ 



(20) 



Thus, in the limit as fig — > oo the generalized potential takes the simple form 



U{q 2 ) = U L {q 2 ) + C with C = T{-fi 2 M ) - U L (-fi 2 M ) 



(21) 



at least in the region, q 2 > —fi 2 M , relevant for (12). The potential Ui(q 2 ) is 
the one implied by the long-range force only. It was studied in some detail in 
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Fig. 5. Case study with repulsive short-range force. The l.h.p. and r.h.p. follow with 
repulsive and attractive long-range potential. 

the previous section. The size of the coupling constant g$ can be used to dial 
the value of the scattering amplitude at the matching point. This is nothing 
but the desired renormalization condition which we will be insisting on in the 
following. Owing to this condition one may view our scheme as an analytic 
continuation of the scattering amplitude from the matching point onwards 
into the physical region. 

We consider the case of a repulsive and attractive short-range interaction with 
fj,s — 12 [Ml for defmiteness. At the scale of the long-range potential the short- 
range potential can be viewed as pointlike. We set 



|^| = 0.95^ 5 ~0.80^, (22) 

m 

so that the interaction is rather strong, but a bound state is not yet formed 




Fig. 6. Case study with attractive short-range force. The l.h.p. and r.h.p. follow 
with repulsive and attractive long-range potential. 
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Fig. 7. Binding energy for 
the system with two attractive 
Yukawa interactions as a func- 
tion of g^. 
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for the attractive case. For the long-range part of the potential we take the 
coupling constants g — ±^f. The influence of the left-hand cut implied by 
the short-range potential can be further weakened by a correlation of the 
parameter A and the Yukawa mass fis- For 



A= 2 fis 



(23) 



a smooth spectral weight at the leading branch point with q 2 = —fi 2 s / 4 arises. 
Any different choice leads to a logarithmic divergence at that branch point. 



In the following we use (23) as it implies a significantly more rapid realization 



of (21 ) as fis — >■ oo. 



Given the parameter set we compute the phase shift for the four different 
scenarios. In Fig. [5] and Fig. [6] the exact results are confronted with the ap- 
proximation implied by (21). The shaded area in the figures illustrate the 



accuracy of (21). While for the repulsive short-range force the thick and thin 



solid lines are almost indistinguishable, we observe small discrepancies for the 
attractive short-range force. The thin lines correspond to the approximate re- 



sults with (21). The additional dashed and dashed-dotted lines follow with 
the long-range potential Ui(q 2 ) in (21) truncated according to (17). In all 



cases already the second order provides a reasonable approximation, where 
the convergence is more rapid for the repulsive short-range force. 

We further scrutinize the quality of our approximation scheme. Fig. 7 shows 
the dependence of the binding energy on g L < — 0.11g cL of the bound state 
that is formed for the attractive short-range potential and a sufficiently strong 
long-range force. Shown is the exact binding energy together with the approx- 



imations defined by (21) and (17). The solid lines are as in the two previous 



figures. The determination of the binding energy is quite accurate already at 
zeroth order in the conformal mapping for a wide range of couplings constants. 



We summarize that the constraints set by causality and unitarity can be used 
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to arrive at a quite effective expansion scheme that is suitable for applications 
in effective field theories pQ. 
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